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We present a measurement of the branching fraction of the decay B° ~ * a* (1260) 
with af (1260) — ► ■k^tt ± h ± . The data sample corresponds to 218 x 10 6 BB pairs pro- 
duced in e + e~ annihilation through the T(4S) resonance. We measure the branching frac- 
tion B(B° -> of (1260) -k^)B{ci± (1260) -» 7r T 7r ± 7r ± ) = (16.6 ± 1.9 ± 1.5) x 10~ 6 , where the first error 
quoted is statistical and the second is systematic. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The rare decay B° -> af (1260) ir T is expected to be 
dominated by b — > uud contributions. For the branching 
fraction of this decay mode an upper limit of 49 x 10~ 5 
at the 90% C.L. has been set by CLEO Q. Bauer et 
al. have predicted a branching fraction of the decay 
B° -> a~(1260) 7T+ of 38 x 10~ 6 within the framework of 
the factorisation model and assuming |T4fe/^cb| = 0-08 @. 
The study of this decay mode is complicated by the large 
discrepancies between the parameters of the ai(1260) me- 
son obtained from analyses involving hadronic interac- 
tions 0] and r decays 0. The decay B° -> of (1260) tt^, 
in addition to the decays B° — ► 7r + 7r~, B° — > /9 ± 7r =F , and 
B° — > p + p~ , can be used to give a new measurement of 
the Cabibbo-Kobayashi-Maskawa angle a of the Unitar- 
ity Triangle 0- 

We present a measurement of the branching fraction of 
the decay B° -> af (1260) tt^ with of (1260) -> 7rT 7 r ± 7r ± . 
The ai(1260) — > 37r decay proceeds mainly through the 
intermediate states (7r7r) p 7r and (tvtv)^ |6j. No attempt 
is made to separate the contributions of the dominant P- 
wave (7T7r)p and the S-wave (7r7r) CT in the channel 7r + 7r~. 
Only a systematic uncertainty is estimated due to the 
difference in the selection efficiency. Possible background 
contributions from B° decays to B° — > af (1320) ir T and 
B° — > 7r ± (1300) tt t are investigated. 

The data were collected with the BaBar detector Q at 
the PEP-II asymmetric e + e _ collider 8]. An integrated 
luminosity of 198 fb" 1 , corresponding to 218 million 
BB pairs, was recorded at the T(4S I ) resonance ("on- 
resonance", center-of-mass energy yfs = 10.58 GeV). 
An additional 15 fb _1 were taken about 40 MeV below 
this energy ( "off-resonance" ) for the study of continuum 
background in which a light or charm quark pair is pro- 
duced instead of an Y(4S). 

Charged particles are detected and their momenta 
measured by the combination of a silicon vertex tracker, 
consisting of five layers of double-sided silicon detectors, 
and a 40-layer central drift chamber, both operating in 
the 1.5-T magnetic field of a superconducting solenoid. 
The tracking system covers 92% of the solid angle in the 
center-of-mass frame. 

Charged-particle identification (PID) is provided by 
the average energy loss (dE/dx) in the tracking devices 
and by an internally reflecting ring-imaging Chcrcnkov 
detector (DIRC) covering the central region. A K/tt 
separation of better than four standard deviations (a) 
is achieved for momenta below 3 GeV/c, decreasing to 
2.5 er at the highest momenta in the B decay final states. 



Monte Carlo (MC) simulations of the signal decay 
modes, continuum, BB backgrounds and detector re- 
sponse [9j are used to establish the event selection crite- 
ria. The MC signal events are simulated as B° decays to 
ai(1260)7r with ai — > pn. For the ai(1260) meson param- 
eters we take the mass too = 1230 MeV/c 2 and Tq = 400 
MeV/c 2 0,E3 ■ 

We reconstruct the decay af (1260) — > 7r =F 7r ± 7r ± with 
the following requirement on the invariant mass: 0.83 < 
m ai(1260) < GeV/c 2 - The intermediate dipion state 
is reconstructed with an invariant mass between 0.51 and 
1.1 GeV/c 2 . We impose several PID requirements to en- 
sure the identity of the signal pions. For the bachelor 
charged track we require an associated DIRC Chcrcnkov 
angle between —2 a and +5 a from the expected value 
for a pion. With this requirement all but 1.4% of any 
background from ai(1260)-ftT is removed. 

A B meson candidate is characterized kine- 
matically by the energy-substituted mass toes — 
\J (s/2 + p • pb) 2 /Eq — p' 2 B and energy difference 
AE = Eg — ^y/s, where the subscripts and B refer 
to the initial T(4S) and to the B candidate in the 
lab-frame, respectively, and the asterisk denotes the 
T(4S*) frame. The resolutions in toes and in AE are 
about 3.0 MeV/c 2 and 20 MeV respectively. We require 
\AE\ < 0.2 GeV and 5.25 < to es < 5.29 GeV/c 2 . To 
reduce fake B meson candidates we require a B vertex 
X 2 probability > 0.01. The cosine of the angle between 
the direction of the tt meson from ai(1260) — > pn with 
respect to the flight direction of the B in the ai(1260) 
meson rest frame is required to be between —0.85 
and 0.85 to suppress combinatorial background. The 
distribution of this variable is flat for signal and peaks 
near ±1 for this background. 

To reject continuum background, we use the angle Qt 
between the thrust axis of the B candidate and that of 
the rest of the tracks and neutral clusters in the event, 
calculated in the center-of-mass frame. The distribu- 
tion of cos 9t is sharply peaked near ±1 for combinations 
drawn from jet-like qq pairs and is nearly uniform for the 
isotropic B meson decays; we require | cos 6t\ < 0.65. 
The remaining continuum background is modeled from 
off-resonance data. We use MC simulations of B°B° and 
B + B~ decays to look for BB backgrounds, which can 
come from both charmless and charm decays. We find 
that the decay mode B° — * D~ir + , with D~ — > K + ir~ir~ 
or D~ — > Kgir^ , are the dominant BB backgrounds 
to ultimate final states different than the signal. The 
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decay modes B° -> af (1320)ttF and B° -> ^ (1300) ^ 
have the same final daughters as the signal. We sup- 
press these with the angular variable A, defined as the 
cosine of the angle between the normal to the plane of 
the 37r resonance and the flight direction of the bachelor 
pion evaluated in the Sir resonance rest frame. Since the 
£H(1260), a 2 (1320) and vr(1300) have spins of 1, 2 and 
respectively, the distributions of the variable A for these 
three resonances differ. We require \A\ < 0.62. 

We use an unbinned, multivariate maximum- 
likelihood fit to extract the yields of B° — > af (1260) ir^, 
B° -> af (1320) ttT and B° -> w ± (1300) tt^. The likeli- 
hood function incorporates five variables. As mentioned 
above, we describe the B decay kinematics with two 
variables: AE and i«es- We also include the invariant 
mass of the 37r system, a Fisher discriminant T 1 and 
the variable A (though the later provides little discrim- 
ination after the requirement mentioned above). The 
Fisher discriminant combines four variables: the angles 
with respect to the beam axis, in the T(4S*) frame, of 
the B momentum and B thrust axis, and the zeroth and 
second angular moments L0.2 of the energy flow around 
the B thrust axis. The moments are defined by 

Lj = x Icos^l 3 , (1) 

i 

where 8i is the angle with respect to the B thrust axis 
of track or neutral cluster i, pi is its momentum, and 
the sum excludes tracks and clusters used to build the B 
candidate. 

We have on average 1.4 candidates per event and we 
select the B candidate with the smallest \ 2 formed from 
the p mass. The efficiency of the best candidate algo- 
rithm is 94%. 

Since the correlation between the observables in the se- 
lected data and in MC signal events is small, we take the 
probability density function (PDF) for each event to be 
a product of the PDFs for the separate observables. The 
product PDF for event i and hypothesis j, where j can be 
signal of (1260) tt^, af (1320) 7r=F and ?rf (1300) ?rF back- 
grounds, continuum background or BB background (2 
types), is given by: 

V) = Vj (m ES ) • Vi (AE) ■ Vj (T) ■ V 3 (m ai )-V j {A). (2) 

The probability that inside the signal event the primary 
pion from the B candidate is confused with a pion from 
the di(1260) is negligible because of the high momentum 
of the primary pion in T(45) frame. There is the possi- 
bility that a track from a af (1260) ir T , af (1320) ir T and 
7r ± (1300) 7r T event is exchanged with a track from the 
rest of the event. These so-called self cross feed (SCF) 
events are considered as background events. The likeli- 



hood function for the event i is defined as 

C = J2 {"M + n s k CF V s k CF ) + n qq V qq (3) 

k=l 

+ n BB\P BB\ + n BB2'BB2 > 

where n k and nf (k = 1, 3) are the signal and SCF 
yields for af (1260) tt^, af (1320) tt^, and 71^(1300) ttF, 
respectively, n qq is the number of continuum background 
events, n BB1 is the number of BB background events 
D~tt + with D~ — » K + tt~t:~ and n BB2 is the number 
of BB background events D~tt + with D~ — > Kgir~. 
Vk is the PDF for correctly reconstructed MC signal 
events;?^ CF is the PDF for SCF events, V qq is the PDF 
for continuum background events, and V b q x and V BB2 
are the PDFs for the two types of BB backgrounds, all 
evaluated with the observables of the ith event. 

We write the extended likelihood function for all events 
as : 

N 

£ = exp(-^n J )n^' ; (4) 

j i 

where rij is the number of events of hypothesis j found by 
the fitter, and N is the number of events in the sample. 
The first factor takes into account the Poisson fluctua- 
tions in the total number of events. 

We determine the PDFs for signal and BB back- 
grounds from MC distributions in each observable. For 
the continuum background we establish the functional 
forms and initial parameter values of the PDFs with off- 
resonance data. 

The PDF of the invariant mass of the ai(1260) meson 
in signal events is parameterized as a relativistic Breit- 
Wigner line-shape with a mass dependent width which 
takes into account the effect of the mass-dependent p 
width [h]. The PDFs of the invariant masses of the 
a2(1320) and 7r(1300) mesons are parameterized by triple 
Gaussian functions. The tties and AE distributions 
for signal are parameterized as double Gaussian func- 
tions. The AE distribution for continuum background 
is parameterized by a linear function. The combinato- 
rial background in mES is described by a phase-space- 
motivated empirical function [l^. We model the Fisher 
distribution T using a Gaussian function with different 
widths above and below the mean. The A distributions 
are modeled using polynomials. 

In the fit there are fourteen free parameters: six yields, 
the signal ai(1260) mass and width, and six parameters 
affecting the shape of the combinatorial background. Ta- 
ble I lists the results of the final fits. Fitted values of 
ai(1260) mass and width have statistical errors only. 

Equal production rates to B°B° and B + B~ pairs 
are assumed. We find no evidence of the decay 
B° -> 7^(1300) ttF, and therefore we have not included 
this component in the final fit. The yield of the decay 
B° -> af (1320) ?rF is 8.3 ± 23.6 events. 
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TABLE I: Signal yield, detection efficiency (e), statistical sig- 
nificance (with systematic uncertainties), branching fraction, 
and the mass and width of the ai(1260) meson. 



Fit quantity 


af (1260) tt t 


Signal yield 


421 ± 48 


e (%) 


11.7 


Stat. sign, (a) 


9.2 


6(xHT 6 ) 


16.6 ± 1.9 ± 1.5 


m(ai(1260)) 
r(ai(1260)) 


1229 ± 21 MeV/c J 
393 ± 62 MeV/c 2 




-0.2 -0.1 0.1 0.2 5.25 5.26 5.27 5.28 5.29 
AE (GeV) m ES (GeV / c 2 ) 




m a (GeV / c 2 ) Fisher 

FIG. 1: Projections of a) AE, b) tties, c) m ai , and d) T. 
Points represent on-resonance data, dotted lines the contin- 
uum and BB backgrounds, and solid lines the full fit func- 
tion. These plots are made with a cut on the signal likelihood 
which includes about 40% of the signal. 

We find a signal yield bias of +3.8% by generating 
and fitting MC simulated samples containing signal and 
background populations expected for data. We find that 
— In L max from the on-resonance data lies well within the 
distribution of — lnL max from these simulated samples. 
The signal reconstruction efficiency is obtained from the 
fraction of signal MC events passing the selection criteria, 
adjusted for the bias in the likelihood fit. The statistical 
significance is taken as the square root of the difference 
between the value of —2 In L for zero signal and the value 
at its minimum. 

In Fig. H we show the AE, tties, Tn ai , and T projec- 
tions made by selecting events with a signal likelihood 
(computed without the variable shown in the figure) ex- 
ceeding a threshold that optimizes the expected sensitiv- 
ity. The enhancement at 1.7 GeV/c 2 in Fig. [He) comes 
from D~it + background. 

In Fig. [21 we show the distribution of the ratio of the 
likelihood for signal events L(Sg) and the sum of likeli- 
hoods for signal and all types of background [L(Sg) + 
L(Bg)] for on-resonance data and for Monte Carlo events 
generated from PDFs. We see good agreement between 
the model and the data. By construction the background 
is concentrated near zero, while the signal appears as an 



excess of events near one. 




1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
L(Sg)/[L(S 8 ) + L(Bg)] 



FIG. 2: Likelihood ratio L(Sg)/[L(Sg) + L(Bg)]. Points rep- 
resent the data, the solid histogram is from Monte Carlo sam- 
ples of background plus signal, with the background compo- 
nent shaded. 

Most of the systematic errors on the signal yield that 
arise from uncertainties in the values of the PDF pa- 
rameters have already been incorporated into the overall 
statistical error, since they are floated in the fit. We 
determine the sensitivity to the other parameters of the 
signal and background PDF components by varying these 
within their uncertainties. 

The uncertainty in our knowledge of the efficiency is 
found to be 3.2%. The systematic error on the fit yield is 
6.2%, which is obtained by varing the PDF parameters 
within their uncertainties. We estimate the uncertainty 
in the number of BB pairs to be 1.1 %. The uncertainty 
in the fit bias correction is 1.9%, taken as half of the 
fit bias correction. Published world averages 0| provide 
the B daughter branching fraction uncertainties. The 
systematic errors on 01(1260)^ cross-feed background 
and on SCF are both estimated to be 1.4%. The po- 
tential background contribution from B° decays to p°p°, 
p°ir + ir~ and 47r is estimated assuming the branching 
fractions of 1, 2, and 2 in 10~ 6 respectively [ljj. The 
associated systematic uncertainty is 3.9%. The system- 
atic effect due to differences between data and MC for 
the cos6*t selection is 1.8%. A systematic uncertainty 
of 2.5% is estimated for the difference in reconstruction 
efficiency in the decay modes through the dominant P- 
wave (7T7r) p and the S-wave (tttt)^. The contribution of 
interference between ai (1320) and ai(1260) is negligible. 
In fact, varying the O2(1320)7r background with different 
selection criteria on the angular variable A gives no sig- 
nificant change to the efficiency-corrected signal yield of 
ai(1260)7r. We find also that the systematic effect due 
to different form factors in MC signal simulation is neg- 
ligible. The total systematic error is 9.1 %. 

In conclusion, we have measured the branching frac- 
tion B(B° -> a±(1260) 7rF)2?(af (1260) -> 7r = F 7 r ± 7r ± )= 
(16.6 ± 1.9 ± 1.5) x 10" 6 . Assuming B(af (1260) -> 
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7rT7r ± 7r ± ) is equal to B(af (1260) -> 7r ± 7r°7r°), and that 
B(af (1260) -> (3tt) ± ) is equal to 100% @], we obtain 
B(B° -> of (1260)^) = (33.2 ± 3.8 ± 3.0) x 10~ 6 
The decay mode, observed for the first time, is seen 
with a significance of 9.2 a, which includes systematic 
uncertainties. 
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